Reheating of the universe after inflation with f(4>)R gravity 
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We show that reheating of the universe occurs spontaneously in a broad class of inflation models 
with f(<j))R gravity ((j> is inflaton). The model does not require explicit couplings between (j> and 
bosonic or fermionic matter fields. The couplings arise spontaneously when <j> settles in the vacuum 
expectation value (vev) and oscillates, with coupling constants given by derivatives of f((j>) at the 
vev and the mass of resulting bosonic or fermionic fields. This mechanism allows inflaton quanta to 
decay into any fields which are not conformally invariant in f(<j>)R gravity theories. 



Inflation is an indispensable building-block of the stan- 
dard model of cosmology [l], , and has passed a number 
of stringent observational tests Any inflation mod- 
els must contain a mechanism by which the universe re- 
heats after inflation [4j. The reheating mechanism re- 
quires detailed knowledge of interactions (e.g., Yukawa 
coupling) between inflaton fields and their decay prod- 
ucts. Since the physics behind inflation is beyond the 
standard model of elementary particles, the precise na- 
ture of inflaton fields is currently undetermined, and the 
coupling between inflaton and matter fields is often put 
in by hand. On the other hand, inflaton and matter fields 
are coupled through gravity. Of course the gravitational 
coupling is suppressed by the Planck mass and thus too 
weak to yield interesting effects [l|; however, we shall 
show that the reheating occurs spontaneously if inflaton 
is coupled to gravity non-minimally, i.e., the gravitational 
action is given not by the Einstein-Hilbert form, R, but 
by f(4>)R, where <fi is an inflaton field. Even if matter 
fields do not interact with <f> directly, they ought to inter- 
act via gravitation whose perturbations lead to Yukawa- 
type interactions. A similar idea was put forward by Q, 
who considered preheating with a non-minimal coupling 
between matter and gravity. Here, we do not consider 
preheating, but focus only on the perturbative reheating 
arising from f((f))R gravity. Therefore, we can calculate 
the resulting reheating temperature after inflation. 

Why study _J{4>)R gravity? There are a number of mo- 
tivations 0, 0]. The strongest motivation comes from 
the fact that almost any candidate theories of funda- 
mental physics which involve compactification of extra 
dimensions yield f{<j))R with the form of /(</>) depending 
on models. One illuminating example would be string 
moduli with f(<f)) cx e~ a<t> . Zee's induced gravity theory 
8] has f(4>) = £(f) 2 , and renormalization in the curved 
spacetime yields other more complicated higher deriva- 
tive terms @. Classic scalar-tensor theories, originally 
motivated by Mach's principle [lOj , also fall into this cat- 
egory. It has been shown that inflation occurs naturally 



in these generalized gravity models [Ill.ll2l. ll3lll4ll5l . llrj| 

and the spectrum of scalar curvature perturbations [171 ] 
as well as of tensor gravity wave perturbations [l8| can be 
affected by the presence of f(<fi)R, thereby allowing us to 
constrain f((f>) from the cosmological data. We show that 
limits on the reheating temperature (e.g., gravitino prob- 
lem) provide additional, totally independent constraints 
on f{4>). This argument opens up new possibilities that 
one can constrain a broad class of inflation models from 
the reheating temperature. Reheating in Starobinsky's 
R 2 inflation has been considered by [la, [l!| . 

Interactions between <f> and matter stem from a mixing 
between metric and scalar field perturbations through 
f((f>)R 6}- It is this "gravitational decay channel" from 
(j) to matter that makes the universe reheat after inflation. 
We realize this in a simple Lagrangian given by 



+ An, (1) 



where C m is the matter Lagrangian. Note that there is 
no explicit coupling between 4> (inflaton) and C m . To 
guarantee the ordinary Einstein gravity at low energy, 
we must have f(v) — Mpj, where v is the vacuum ex- 
pectation value (vev) of <fi at the end of inflation and 
Mpi = (87TG)" 1 / 2 = 2.436 x 10 18 GeV is the reduced 
Planck mass. This Lagrangian satisfies the weak equiv- 
alence principle: (f> does not couple to matter directly, 
v = 0- This, however, does not imply that in- 
flaton quanta, i.e., fluctuations around the vev, cannot 
decay into matter. The gravitational field equation from 
Eq. H]) is 



f(4>)G, 



T% + d^d v cj> - -g^{dct>f - g^Vfr) 



(2) 
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Before reheating completes the energy density of the uni- 
verse was dominated by (j>; thus, we treat T™ v as perturba- 
tions. As usual we decompose t/ M „ into the background, 
g^v, and perturbations, h^, as g^ = + h^. The 
first-order perturbation in G^ v is given by 
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where W v = g fJ,x h\ u and h = g^h^. 

Reheating occurs at the potential minimum of <f> where 
<j> oscillates about v. We thus expand <p as <f> = v + 
a, where a represents inflaton quanta which decay into 
matter fields. Treating a as perturbations, one obtains 
the linearized field equation: 



/(»)&„□ " V„V„)<7 = T- (4) 



where f'(y) means df /d4>\tf, =v . This equation contains 
both nhpv and Dcr, and thus wave modes are mixed up 
together. To diagonalize the wave mode we define a new 
field as [| 



(5) 



where the tilde denotes the operation defined by this 
equation, which essentially corresponds to the infinitesi- 
mal conformal transformation. The inverse operation is 



) -h I- l-l^l- 



(6) 



Using the harmonic gauge (analog of the Lorenz gauge) 
conditions, Va/i a „ = 0, one can show that h pv indeed 
obeys the linearized field equation for the wave mode 
(see, e.g., Eq. (35.64) in [H; Eq. (10.9.4) in [Hj]) 



Hhfiv — 2R il \ pl jh Xp — 2R p ^h p ^ — g pv R\ p h Xp + Rh^ 

2 rpm 

~ T\/f2 1 IW 



(7) 



One may also impose the traceless condition on h pu to 
make sure that h pv describes tensor (spin two) gravity 
waves, h — —h — 4/'(w)cr/M|j = 0, which relates a trace 
part of the original metric perturbations to a as h = 
—Af'(v)a/Mp V In this sense a describes the "scalar (spin 
zero) gravity waves" , which are common in scalar-tensor 
theories of gravity. For generality we keep h explicitly 
throughout the paper. 

We consider both fermionic, ip, and bosonic, x, fields 
as matter: C m = + C x + Ant- First, let us write C$ 
as 



^gip [Ip + m-4,] ip, 



(8) 



where Ip = e^jaD^, is a tetrad (vierbein) field, and 
D „ = — is the covariant derivative for spinor fields 



l2l| . T p is a spinor connection and S a/3 are generators 
of the Lorentz group given by T AI (x) = — ^ al3 e a x S7 p e/s\ 
and£ Q/3 = -E^ Q = i[7 a ,7 /3 ], respectively. Here a, (3,... 
denote Lorentz indices while y,, v, . . . denote general co- 
ordinate indices. Note that we have not antisymmetrized 
the first term in Eq. |5]) to make the expression simpler. 



The antisymmetrized term yields the same result. Note 
also that we shall ignore the existence of torsion. 

We expand the tetrad and metric into the background 
and perturbations as e^ 01 ~ e^" — e Xa h p, \/2, and \f—g = 
e + Se ~ e (1 + g p -"h fllJ /2), respectively. The Lagrangian 
becomes (with Eq. ©) 

C^p ~ -eip [e^jaDp + m^] ip 

+ \e^[e lJa la {h p v + -W^vJDp + hm^ip 

+eg^aipip. (9) 

The terms that are proportional to h may be set to vanish 
by gauge transformation. Here we have used the back- 
ground Dirac equation, e^jaD^tp = — m^ip, to obtain 
the last term and ignored the second order terms as well 
as thermal mass induced by quantum corrections in ther- 
mal bath [HHH. The last term in Eq. © is a Yukawa 
interaction term with a coupling constant given by 



2M|i ' 



(10) 



Therefore, a can decay into ip and ip. Note that the 
Yukawa interaction vanishes when f'(v) = or = 0, 
which is consistent with previous work (l3| : as massless 
fcrmions are conformally invariant, no peculiar effects 
can be caused by f((p)R gravity at the tree level. Quan- 
tum corrections such as a conformal, or trace, anomaly 
might open additional decay channels. 

Next, let us consider the bosonic matter, x, given by 



(ii) 



which may be expanded as 



-s^x^x 



U(x) + \h^dp X duX 



f'(v) 
2Ml x 



,(12) 



g(l-h/2-2af'(v)/M^) 



where we have used \J —g — 
and gi*" ~ g^ - + g^h/2 + g^af'(v)/M^. Again, 
the terms that are proportional to h may be set to van- 
ish by gauge transformation. For simplicity we assume 
that x is a massive field with self-interaction, U(x) — 
m^x 2 / 2 + Ax 4 /4. The last term in Eq. (12]) then yields 
the following interactions: y/ = g[af , (v)/2Mp^(2myX 2 + 
^X 4 + g^dpxdvX)- We rewrite these interactions as 



-9- 



Kx 2 +9^^Xd v x)) » V^5x^X 2 ,(13) 



where we have defined a coupling constant, 
/'(«) 



9x 



2M^ 
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(14) 
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which will give a rate of a decaying into two vs- To derive 
Eq. US]) we used g^B^X = -xU'(x)+9^ p(x&x), 
where V M is the covariant derivative on the background 
metric, and estimated as crV M (x<9^x) ~ (JW ctX 2 /2- This 
estimation is valid in the zero temperature limit (see [24j 
for the high temperature limit). While the first term 



in Eq. (|14[) is absent when m x = 0, a can still decay 
through the second term as minimally coupled massless 
scalar fields are not conformally invariant. 

What is the physical interpretation of these couplings? 
Through the field mixing that can be seen in Eq. {5{ the 
inflaton quanta are coupled to the trace (spin zero) part 
of h^ v , which describes the scalar gravity waves. The 
scalar waves are then coupled to the matter fields. If 
inflaton quanta are at least twice as heavy as ^ or x, 
they decay into tf>tp or two xs. 

In order to compute decay rates, one must canoni- 
cally normalize fields. 1 The kinetic term of a is given 
by ££ n A/=5 = -U d(J ) 2 [l + |(/»/Mpi) 2 ] to the 
leading order, 



and thus a = 



a 



'l + f(/»/Af P1 ) 2 
is the canonically normalized inflaton field. There- 
fore, coupling constants should be replaced accord- 
ingly by 9i , = 94 , [l + f(/»/M P1 ) 2 ]~ 1/2 and g x = 

9 X [i+§(/»/M P1 ) 2 r i/2 ._ 

The decay rates of a — > ipip from the last term in Eq. @ 
and a — > yy from Eq. (fT3|) , r to t 
given by [25| 
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m 2 



3/2 



tanh(^)c^,(15) 



1/2 



< "lh('^)c x , (16) 



where m 2 = V"(v) and we have included suppression of 
decay rates due to thermal mass, and C x , which are 
unity in the low temperature limit (T -C m a ) in which 
thermal mass is small, but can be quite small otherwise. 
The exact form of C depends on how decay products 
are thermalized. Yokoyama has calculated C x from ther- 
malization due to AU = Ax 4 /4: for T ^> m a /V\ he 
finds C x = X/(8tt 2 ) for decay via ax 2 interaction [23| 
and C x = A/(10247r 2 ) for <j(dx) 2 (which corresponds to 
the last term in Eq. (|12p : if this term dominates g x is 
given by g x = f'(v)XT 2 /(2Mp l ) because of dominance of 
thermal mass, \/XT, over the intrinsic mass of x) [HI; 
while Cif) has not been calculated explicitly yet. 

As usual we use the condition, 3H(t*) = T to t, to 
define the reheating time, t* , at which radiation be- 
gins to dominate the energy density of the universe. 
From the Friedmann equation, H 2 = p/^Mpj), we get 



1 We are grateful to N. Kaloper for pointing out importance of the 
normalization of a fields. 



p(t*) = r 2 ot M 2 !/3 = foff*( r ''h) r rh; where the last equ- 
ity holds if all the decay products interact with each other 
rapidly enough to achieve thermodynamical equilibrium 
[H. if some fraction of inflaton energy density is con- 
verted into the species that never interact with the vis- 
ible sector, the reheating temperature may be lowered. 
g*(T v h) is the effective number of relativistic degrees of 
freedom at the reheating time. Since r to t 
the reheating temperature is bounded from below 



yr tot A/ P i ( g*(T Th ) 

(IOtt 2 ) 1 ^ 



100 



-1/4 



> 



o , 5 2 \ / m a Mpi ( g* (T rh ) 



^ ml 4tt(40)V« 



100 



-1/4 



(17) 



where we have used Eq. ()15|) and ()16|) assuming for sim- 
plicity that decay products are much lighter than a (m^, , 
m x <C m a ), the reheating temperature is much smaller 
than m a (T <C m a ), and therefore suppression due to 
thermal mass is unimportant (C^, = 1 = C x ). General- 
ization to the other limits is straightforward. One may 
reverse the argument and set a limit on f'(v) as 



< 87r(40) 1/4 T rh 



Mm 



PI 



3/2 



100 



1/4 



a«) 



which provides non-trivial constraints on inflation 
models with f(4>)R gravity. Let us consider a 
non-minimal coupling model, f((f>) = M 2 + £(f> 2 , 
for instance (M 2 is given by the condition 
f(v) = AfpJ. We obtain |f|(l + &^ 2 v 2 /M 2 ^ 1 / 2 < 
4^(40) 1 / 4 (T rh /«)(M P1 /m ff ) 3 / 2 ( ff 4T rh )/100) 1 / 4 , which 
provides a new constraint on £, independent of the 
existing constraints from homogeneity and isotropy fl2| 
and curvature as well as tensor perturbations (l8| . 

Finally, let us show that one can obtain the same re- 
sults by performing the following conformal transforma- 
tion, g^ v = VL 2 g^, with 



n 2 = 



M 2 X 



= 1 



f'(v)a , f"{v)a 2 



M 2 , 



2Af 2 ! 



+ 0{a 3 ). (19) 



Hereafter we put hats on variables in the Einstein frame 
(E frame). We use Maeda's transformation formula [26[ 
to obtain the Lagrangian in the E frame, 



£ = V-g 



M 2 - 1 ~ - A - 

-^lR--g^d^d^-V( ( f > ) 



+ C m , (20) 



where R is calculated from g llVl a transformed poten- 
tial is given by V(<fi) = r2 _4 V(</>), and a new scalar 
field, <p, is defined such that it has the canonical ki- 
netic term. It is related to the original field by d(f)/d<f> = 

Mny/l/ftt) + | (/W/(0)) 2 . As 4> is minimally cou- 

pled to R, there is no mixing between <? M „ and 4>. 



4 



A transformed Lagrangian for ip in the E frame is 

C 4 , = -ei]j[$ - $)/2 + r'm^, (21) 

where yb = ft" 3 / 2 V>, 4> = fr 3/2 ?/S, e» a = {l~ 1 e fJ " a , and 
e = f2 4 e. While no Yukawa interaction arises from the 
kinetic term, the mass term yields a Yukawa interaction: 



e— W 



2M|i 



(22) 



Thus, Yukawa coupling constants agree precisely. 
A transformed Lagrangian for x m the E frame is 



gg^xDvX 



9U(x), (23) 



where x = 1 x and we have followed the procedure in 
6] to define a "covariant derivative" and potential in the 
E frame as = + d^lnQ), and L>(x) = ^~ 4 t/(x), 
respectively. In the case of U(x) = m 2 x 2 /2 + Ax 4 /4, the 
second term of Eq. (|2"3")) becomes 



2 ™^ 2 + 4 X 4 2M2 



^X 



pi 



(24) 



to the linear order in a; the last term agrees with the 
first term in Eq. (I13|) . The covariant derivative yields 
a coupling, —^/—gxg^i^Xii^ hi£l). After integration 
by parts and lnf2 ~ [/'(w)/(2Afp 1 )]cr this term yields 



coupling constants are proportional to the higher order 
derivatives such as f"(v), f"'(v), etc., actually domi- 
nate if /((/>) also has a minimum at the vev, f'(v) 
0. We must canonically normalize quanta as a — 

cry 1 + |(/'(u)/Mpi) 2 when calculating their interaction 

rates. While we have ignored a parametric resonance 
(preheating) [27j entirely, it would be interesting to study 
how preheating might occur in the present context. 

In summary, we have presented a natural mechanism 
for reheating of the universe after inflation without in- 
troducing any explicit couplings between inflaton and 
fermionic or bosonic matter fields. This mechanism al- 
lows inflaton quanta to decay into any fields which are 
present at the end of inflation and are not conformally in- 
variant, when inflaton settles in the vacuum expectation 
value and oscillates. Reheating therefore occurs sponta- 
neously in any theories of f{<j))R gravity. 

We have calculated the reheating temperature from 
this mechanism. We argue that one must always check 
that the reheating temperature in any f((f>)R inflation 
models is reasonable, e.g., the reheating temperature 
does not exceed the critical temperature above which 
too many gravitinos would be produced thermally. This 
mechanism might also allow <f> to decay into graviti- 
nos. How our mechanism is related to an inflaton decay 
through supergravity effects [28| merits further investi- 
gation. Both of these effects should provide non-trivial 
constraints on f(<fr) from the cosmological data. 



(M-vdviV^xg^d^x) -- 

(25) 

which agrees with the second term in Eq. (|13p precisely. 

One can calculate the other interaction terms such as 
cr 2 x 2 , o^vbvb, tj2 (^x) 2 J etc, with known coupling con- 
stants in the E frame easily. These interactions, whose 
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